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Abstract 

Dynamical reaction theory is presented for synthesis of superheavy elements. Characteristic features 
of formation and surviving are discussed, which combinedly determine final residue cross sections of 
superheavy elements. Preliminary results on Z=114 arc also given. 



1 Introduction 

Superheavy elements around Z = 114(or 126) and N = 184 have been believed to exist according to 
theoretical prediction of stability given by the shell correction energy in addition to average nuclear 
binding energy, jlj This means that heavy atomic nuclei with fissility parameter x > 1 could be stabilized 
against fission by a huge barrier which is resulted in by the additional binding of the shclll correction 
energy around the spherical shape. In other words, if superheavy compound nuclei(C.N.) are formed in 
such high excitation that closed shell structure is mostly destroyed, they have no barrier against fission 
and thus are inferred to decay very quickly. Therefore, the point is how to reach the ground state of the 
superheavy nuclei, or how to make a soft-landing at them. In order to minimize fission decays of C.N. 
or maximize their survival probabilities, so-called cold fusion reactions have been used, which succeeded 
in synthesizing SHEs up to Z = 112. |2) They have the merit of large survival probabilities, but suffer 
from the demerit of small formation probabilities because of the sub-barrier fusion. On the other hand, 
so-called hot (warm) fusion reactions have the merit of expected large formation probabilities and the 
demerit of small survival probabilities due to relatively high excitation of C.N. formed. An optimum 
condition for large residue cross sections of SHEs is a balance or a compromise between formation and 
survival probabilities as a function of incident energy or excitation energy of C.N. formed. || Therefore, 
the whole reaction process has to be described, from the encounter of incident ions to the formation of 
compound nuclei, and then to fission decays with residues left of a small probability. An optimum path 
has to be searched for over all possible incident channels including secondary beams becoming available 
sooner or later. Whether the reaction is "cold" or "hot" fusion is automatically determined by choosing 
combinations of projectiles and targets. The theoretical framework, thus, has to accommodate both 
of them. As is described in section the strong dissipation of energies of nuclear collective motions, 
which is well recognized by the fusion hindrance^ and the long fission life,|| has to be taken into 
account in formation process as well as in survival process. 



2 Theoretical framework 



In reactions of massive systems, it is not clear that the compound nucleus theory can apply, i.e., that 
the formation and the decay of compound nuclei are independent. But for simplicity, we assume that 
is valid, at least approximately. Then, the residue cross section is given by the following formula as 
usual, 

°she = ^E( 2J + !) ' P for(Ec.m.) ■ P s ' urv (E*) (1) 
, where P J t or and Pf url , denote the formation probability and the survival probability, respectively. 



2.1 Formation probability 

Pf or is not equal to a simple transmission coefficient Tj(Ec.m.) of an optical potential, nor to a barrier 
penetration factor Pj(Ec.m.) of the combined potential of Coulomb repulsion and nuclear attraction. 
As is well known, [0 there is the fusion hindrance in massive systems, roughly those with Z1Z2 > 
1,800, which could be interpreted by the overcome of so-called conditional saddle under the strong 
dissipation. The necessary additional energy is so-called extra-push energy. Therefore, Pj or should 
take into account both the usual barrier penetration and the dissipative dynamics over the saddle. The 
former is approximated by the penetration factor for the parabolic barrier, while the latter is treated 
by Langevin equation, concerning the collective shape degrees of freedom of the compound systems 
starting with the contact configuration of the projectile and the target. In one-dimension, Langevin 
eq. is given as 



d 2 q dV dq 

The friction coefficient 7 is taken to be equal to that of so-called one-body wall-and-window 
formula. Q The last time-dependent term R(t) in Eq.(jlj) is a random force associated with the fric- 
tion force. The last two terms are a phenomenological description of effects of the nucleonic degrees of 
freedom in excitation, considered as a heat both. Thus, we assume the fluctuation-dissipation theorem, 



< R(t) ■ R(t') >= 7 • T • S(t - t') (3) 

, where < > denotes an average over all the possible realizations of assumed gaussian noise of R and 
T does the temperature of the heat bath, i.e., the compound nucleus. The delta function in the r.h.s. 
is due to the Markovian assumption. Equivalently, one can use Kramers eq. 

df(q,p,t) d p d dV d d 

—dT- - { -Y q m + d p -^ + dp-^- p + ^T-)}f{ q ,p, t) (4) 

, where p denotes the conjugate momentum to q, and j3 is the reduced friction coefficient 7/m. 
Eq.(^) describes a time-evolutionl of distribution function in the phase space, while Eq. (|2j) does tra- 
jectories of Brownian particle. Actually, Eq.(Q) was proposed by Kramers Q in order to interpret the 
fission width proposed by Bohr and Wheeler, from the dynamical view-point. 
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Figure 1: 
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Figure 2: 



In realistic situations, processes are in many dimensions including mass asymmetry degree of free- 
dom etc. in addition to the elongation or the separation between fragments. An important case that 
we will discuss below is that the incident channel is with Zi-Z2 > 1,800 and the compound nucleus is 
with Z = 114. In Fig. [I], an example of contour maps of energy surfaces is given for the symmetric 
mass partition (ZiZi> 3,000). The axes are the elongation (Zo), i.e., the relative distance between 
fragments and the deformation where fragment deformations are assumed to be proportional to their 
masses. It is seen that the contact point of two spheres is located below the conditional saddle point. 
The system has to climb up, governed by the fluctuation-dissipation dynamics (the multi-dimentional 
version of Eq.(|^))in order to reach the spherical configuration of the total system beyond the saddle 
point. Apparently, an additional incident energy is required, which corresponds to the extra-push 
energy. On the other hand, if we think about a very mass-asymmetric incident channel, such as 48 Ca 
+ 244 Pu (7jiZj2— 1,800), we may not suffer from the extra energy, because the contact point of two 
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incident nuclei can be located inside the conditional saddle point. Actually, the mass asymmetry 
parameter a = 0.66 case is shown in Fig. ^], where the contact point is slightly inside and then the 
composite system evolves mainly over the flat region around the spherical configuration after touching 
each other. This aspect is surely favourable for formation probability. Anyhow, Langevin trajectories 

for- 



are calculated on the energy surfaces in order to obtain Yi 



2.2 Survival probability 

The survival probability P stirt , is usually given as follows, 

Psu ™ (ir) = r n (E*) + T f (E*) (5) 

, assuming that fission and neutron emission are dominant decay modes. The neutron emission 
width r„ can be calculated by Weiskopf formula, while the fission width by Bohr- Wheeler one or 
Kramers' stationary limit. &@ If E* is high enough for more than one neutron emission, then the 
expression at r.h.s. of Eq.(T|) is repeatedly used for sequential E* and factorized to obtain the final 
P SU rv Anyhow, this method is valid only for cases with the fission barrier height being larger than the 
temperature T. In SHE compound nuclei, especially in hot fusion reactions, there might be no fission 
barrier, because the barrier is given by the shell correction energy around the spherical configuration 
or somewhere in deformed region which is expected to melt in high excitation. The potential landscape 
changes depending on the temperature T, i.e., on neutron emission. Then, a dynamical treatment is 
called for. In other words, final residue cross sections of SHE are to be determined by a competition 
between time scale of fission and that of restoration of the shell correction energy due to cooling 
by neutron evaporation. Therefore, a speed of cooling and a temperature dependence of the shell 
correction energy are crucial. The latter is shown in Fig. |^, the energy being normalized by the 
maximum value, i.e., that of T = 0. The calculations are made by the use 
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of single particle levels in a Saxon- Wood finite potential well. The present results are in good agreement 
with the prescription proposed by Ignatyuk. |8) Using this temperature dependence and cooling obtained 
by the statistical code for evaporation, we obtain potential energy curves as a function of time. With 
this time-dependent potential, we calculate fission process with Eq.(||) or (^) or its simplified version 
and obtain a small fraction of probability captured by the pocket made by the restored shell correction 
energy, which gives a final survival probability P sur . 

In Fig. f|, P SMr 's with J=10 are shown for several isotopes of compound nuclei with Z=114, as a 
function of initial excitation energy. We readily see that P sur depends strongly on excitation energy as 
well as on mass number, i.e., neutron number of initial compound nucleus. This is due to differences 
in cooling speed due to different neutron separation energies B„'s. Thus, we have to form a compound 
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nucleus which has as small B„ as possible, i.e., as large neutron number N as possible in order to 
obtain large survival probabilities. In cases with small N such as GSI experiments, higher excitation is 
extremely unfavored, thus "cold" fusion process is used. Note that the results in Fig. |I] are at the time 
2,000 x 10 -21 sec and therefore they should be reduced by a reduction factor due to the last neutron 
emission for the final survival probabilitites. 



3 Preliminary results and discussion 

In brief, using WKB approximation for tunneling and 3-dimensional Langevin equation for subsequent 
dynamical evolutions of the total system, we calculated formation probabilities, while one-dimensional 
Smoluchowski equation which is the approximation of Eq.(|J) was used for fission decay or survival 
probabilities. Evaporations were taken into account in statistical model. Results are shown in Fig. [| 
over possible incident channels leading to Z=114. Existences of the optimum enegies in various channels 
are clearly seen, which resulted in by the compromise of the two factors. Rapid increases in the left 
hand side are due to barrier penetration and dissipative dynamical evolutions of shape degrees of 
freedom, while decreases in the right hand side due to energy dependences of survival probabilities 
against fission. It is worth to notice that the recent Dubna experiment observed an event which would 
be a signature for Z=114 element, at the energy predicted here to be most favourable in the incident 
channel 48 Ca + 244 Pu.[9] 

As for a relation to so-called dinucleus system concept, [^0) it would be helpful to have a look on 
contour map of energy surface in axes of the elongation (Zo), and their mass-asymmetry (a), which 
is shown in Fig. |6|. There is a Businaro-Gallone peak at the upper right corner. Touching points 
with asymmetries a — 0.0 (symmetric case) and a = 0.66 (close to 48 Ca -I- 244 Pu case) are also 
shown, schematically, It should be noted here that the top horizontal line and the left vertical line 
both describe approximately the spherical configuration of the total system, thereby should converge 
to one point in a more realistic representation. The model based on the dinucleus system concept 
gives fusion probabilities over the Businaro-Gallone peak, so it evaluates evolutions along the vertical 
direction in Fig. while the 3-D Langevin calculations take into account evolutions over all possible 
IS^Igij'gi^j^m^a^y including the horizontal direction as well. 
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